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Abstract 

We examine cosmological models with generalized phantom energy (GPE). 
Generalized phantom energy satisfies the supernegative equation of state, but its 
evolution with the scale factor is generally independent, i.e. not determined by its 
equation of state. The requirement of general covariance makes the gravitational 
constant time-dependent. It is found that a large class of distinct GPE models 
with different evolution of generalized phantom energy density and gravitational 
constant, but the same equation of state of GPE have the same evolution of the 
scale factor of the universe in the distant future. The time dependence of the equa- 
tion of state parameter determines whether the universe will end in a de Sitter-like 
phase or diverge in finite time with the accompanying "Big Rip" effect on the 
bound structures. 



Results of recent cosmological observations, such as distant supernovae of type 
la (SNIa) and cosmic microwave background radiation (CMBR) j2j, have dra- 
matically altered our perception of the dynamics and composition of the universe 
and reshaped the landscape of standard cosmology The universe seems to be 
in the phase of accelerated expansion, which started at a relatively small redshift, 
z ~ 1. This acceleration is attributed to a new form of matter, usually referred to 
as dark energy, the nature of which is still not definitely established. Observations 
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indicate that the energy density of the universe is very close to its critical density 
where dark energy presently accounts for approximately 2/3 of the total energy 
density, while the remaining 1/3 comes predominantly from dark matter, another 
unidentified component of the universe. The most prominent and studied candi- 
dates for the title of dark energy are the cosmological constant [U 03 E] (together 
with its dynamical variants, such as renormalization group running cosmological 
constant El El ) , quintessence ^U] and the Chaplygin gas . 

The majority of dark energy models share a common constraint on their equa- 
tion of state (pd and pd represent pressure and energy density of dark energy, 
respectively) 

Pd = wp d , (1) 

where w > — 1. Such a constraint is, however, not justified by the unbiased fits to 
the data of cosmological observations. Moreover, the allowed interval for the pa- 
rameter of the equation of state extends significantly into the region with w < —1. 
The use of observational data on CMBR, large scale structure (LSS), SNIa and 
Hubble parameter measurements from the Hubble Space Telescope (HST) un- 
der the assumption of the redshift independent parameter w give the restriction 
— 1.38 < w < —0.82 at the 95% confidence level ITS]. Therefore, a possible su- 
pernegative equation of state of dark energy deserves due attention. 

A new type of dark energy with the equation of state characterized by w < — 1 
was proposed in and named phantom energy. Phantom energy is considered 
to be separate from other components of the universe and its energy-momentum 
tensor is conserved separately. In such a setting, the equation of state of dark 
energy determines its evolution with the scale factor a. The supernegative nature 
of the equation of state of the phantom energy leads to the growing energy density 
of phantom energy pd ~ a~ s ( 1+w \ for a constant parameter w. The cosmological 
dynamics of the universe with such a phantom energy component possesses many 
interesting features ^1] . The growth of the energy density of phantom energy drives 
the scale factor of the universe to infinity in finite time. The increasing negative 
pressure of phantom energy leads to the unbounding of all bound structures in the 
universe. This dramatic and picturesque scenario of the cosmic doomsday was ap- 
propriately named "Big Rip" . The formulation of microscopic models for phantom 
energy ^3] relies on the machinery developed in quintessence models, namely the 
evolution of the scalar field in a suitably chosen potential. However, the description 
of phantom energy may require an introduction of some nonstandard alterations, 
e.g. the negative kinetic term of the scalar field. Detailed considerations of the 
Lagrangians describing phantom energy show that in some cases the universe with 
phantom energy ends in a "Big Rip" , while in others it asymptotically approaches 
the de Sitter expansion. 

In this paper, we consider models with generalized phantom energy (GPE). 
First, we set up a more general model of the evolution of the universe with phantom 
energy. We assume that there are two components of the universe: the dark 
energy component (which will have the phantom energy characteristics), and the 
"ordinary" matter component with the respective energy densities pd and p m . The 
"ordinary" matter is taken to satisfy the equation of state 

Pm = IPm , (2) 
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where 7 > 0. Furthermore, we assume that the energy-momentum tensor of the 
"ordinary" matter is conserved 

^m\v = . (3) 

The equation given above ensures that the parameter of the equation of state 
governs the evolution of the "ordinary" matter energy density, i.e. 

( a\ -3( 1+ ^) 

Pm = Pm,o — • (4) 

\«0/ 

Dark energy has the equation of state 

Pd = wp d , (5) 

where w generally depends on time explicitly or implicitly, via explicit dependence 
on some other time-dependent quantity, such as the scale factor a. In the case of 
dark energy, we allow the possibility of non-conservation of the energy-momentum 
tensor, i.e. 

T%,*0. (6) 

Thus, the evolution of the dark energy density is not determined by the parameter 
from its equation of state. 

With the properties of the components of the universe defined, we can specify 
the laws of its evolution. We start from the Einstein equation 

CT = StvGT^ , (7) 

where G^ v is the Einstein tensor and T» v = T$ v + is the total energy- 
momentum tensor. The reconciliation of the requirement of the general covari- 
ance of (J7J) and the non-conservation relation @ is possible with the promotion of 
gravitational constant G into a space-time dependent quantity. This change can be 
interpreted as a modification of the dynamics of General Relativity. This additional 
dynamics is effectively described by the introduction of space-time dependence of 
G. We consider the models where G is a function of time only, G = G(t). Mod- 
els with the time-dependent G were extensively studied in the framework of the 
time-dependent cosmological term A(t) JE]- The covariant derivative of ® then 
implies 

(G(t)T^). u = . (8) 
This equation can be rewritten in the form 

d(G(p m + p d )a 3 ) = -G(p m + p d )da 3 . (9) 

Combining the evolution laws (j3J) and Q and introducing w = — 1 + k (where k 
describes the deviation from the parameter of the equation of state inherent to the 
cosmological constant) we arrive at 

G(p m + Pd ) + Gp d + 3KHGp d = 0. (10) 

Here H = a/a is the Hubble parameter, while dots denote time derivatives. Equa- 
tion ()1U|) clearly shows the generality of the model. In the case of the constant G, 
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we recover the standard equation of conservation of Tjf . Equation (|10|) shows that 
the time evolution of G is the result of two competing effects. Namely, for dark 
energy with growing energy density, the second term in 1|1(J|) causes the decrease 
of G, while for negative k, the third term in increases G with time. 

Finally, Friedmann equations for the evolution of the scale factor complete the 
set of evolution equations (jlj) and (fTU|) 

k _ 8vr 
a J 'a 2 3 



7) +72=— G(Pm + Pd), (11) 



(X 47T 

- = — ^-G(p m + p d + 3p m + 3p d ) . (12) 
a 6 

The set of equations (@J), (fTU)l and ((TT|) reveals that we have essentially two 
independent equations for three dynamical quantities G, pd and a (assuming that k 
is the function of these quantities and time). Without a more specific identification 
of the dynamics of G or pd, it is not possible to solve the aforementioned set of 
equations. However, as we show below, with mild assumptions about the evolution 
of dark energy with the scale factor, it is possible to obtain information on the 
future evolution of the universe for general G and pd satisfying the equations given 
above. 

Next, we introduce the concept of generalized phantom energy (GPE). Gener- 
alized phantom energy is the form of dark energy satisfying the equation of state 
(PQ) with the non-conserved energy-momentum tensor ® and the following two 
properties: 

(a) GPE energy density is a non- decreasing function of the scale factor, 

(b) GPE equation of state satisfies k < 0. 

We further examine the future evolution of the universe. In the sufficiently 
distant future we have p m <C pd and p m can be neglected in the evolution equations. 
Equations 1)10(1 and (|llj) thus become 

'd\ 2 k _ 8?r 
. a J a 2 3 



+ - = — Gp d (13) 



and 

j t (G Pd ) + 3 K HGpd = 0. (14) 
Furthermore, from equation (|14j) . we obtain 

^ = -3,> . (15) 
Gp d a 

As the condition — k > is satisfied by assumption (b), we obtain 

Gp d > (Gp d ) • (16) 

Therefore, as Gpd is a growing function in an expanding universe, for large a we can 
disregard the term k/a 2 in equation (|13|) . For the flat universe, this approximation 
is exact, while for the closed or the open universe, this approximation is applicable 
in the sufficienly distant future. 
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Finally, we end up with the following two equations for the dynamics of the 
universe in the distant future: 

H 2 = ^-(Gp d ), (17) 

±(G Pd ) + 3KH(G Pd ) = 0. (18) 
dt 

By combining equations (fT7|) and (|18|). we obtain an equation for the evolution of 
the Hubble parameter H with time 

dH 3 



with the solution 



df +^ 2 = 0, (19) 

H u\ gX^o) (201 

1 + |ff (t ) Jl K(t')dt> ' 

Once we have found the expression for the evolution of the Hubble parameter, it 
is easy to obtain an expression for the evolution of the scale factor a 

a(t) = a(t )exp ( C dt' ^4 j . (21) 

y \Jt l + lH(t )flK(t»)dt») 

General solutions (|2()jl and (|21|) exhibit some interesting features. The evolution 
of the universe in the sufficiently distant future is governed only by the parameter 
of the equation of state of dark energy. The precise form of the growth of p d with 
the scale factor a is irrelevant in this limit. This implies that the entire class of 
models with different functional forms of p d and G, obeying the same equation of 
state, show the same behaviour in the sufficiently distant future. Therefore, we 
can divide all GPE models with the characteristics specified above into classes with 
the same equation of state. 

An important question regarding the fate of the universe is whether, for a 
particular class of generalized phantom energy models, a and H diverge in finite 
time or reach infinite values only in infinite time. For the Hubble parameter H, 
the answer is straightforward. There will be no divergence of H in finite time if 
the denominator of the expression on the right-hand side of Q20J1 remains positive 
for all times. This leads to the condition 

t oil (to) 

As in this case there is no singularity in H (t) in finite time, the scale factor a(t) also 
does not diverge in finite time. In order to have the convergence of the integral 
J^(—K(t'))dt' required in (|22|). the function n(t) has to tend to zero at asymp- 
totically large times. Therefore, for generalized phantom matter which exhibits 
no divergence of H or a in finite time, the parameter of the equation of state 
approaches —1, i.e. generalized phantom energy approaches the time-dependent 
cosmological term. 

In the case when the condition (1221) is not satisfied, the Hubble parameter H 
diverges in finite time t. From Friedmann equations we have 

d = Ha, (23) 
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l-^/tJilV (24) 

These expressions indicate that, when H diverges in finite time t, both a and a 
diverge as well, so the scale factor a cannot remain finite, but diverges in finite 
time t as well. 

From the general expessions (j2U|) and (|21|) . we can obtain evolution laws for 
the conceptually simple, but important case JH] 

K (t) = -k . (25) 

With such a choice for the parameter of the equation of state of generalized phan- 
tom energy, we have the following evolution laws: 



H(t Q ) 

£H(to)Ko(t-to) 



Hit) = ; 3 777^77 ~ > ( 2 6) 
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a{t) = a(t ) (l - |iJ(to)«o(* - to)) 3K ° • (27) 

These solutions clearly show the onset of the divergence in H and a. The universe 
with generalized phantom energy with the constant parameter of the equation of 
state evolves to infinity in finite time. 

Comparison with the case of the "standard" phantom energy j 1 3 1 114j shows 
that, for the same parameter of the equation of state K,(t), the scale factor follows 
the same evolution law. Given the fact that the parameter of the equation of state 
does not determine the scaling with a, and that G is variable in the framework 
of generalized phantom energy, it is by no means obvious that coincidence of this 
sort should exist. However, from the equation (| !()[). we readily see that for the case 
of constant G, we recover the equation of evolution for the "standard" phantom 
energy. As far as the evolution in the sufficiently distant future is concerned, the 
"standard" phantom energy model is just one instance of the class of generalized 
phantom energy models with the same function n(t). 

Given the same evolution properties of the broad class of GPE models with 
the same «(t), it is natural to look at the destiny of bound structures, another 
peculiarity of phantom energy models [12]. The relevant quantity with respect to 
the stability of the bound structures is the analogue of the gravitational potential 
proportional to the quantity G(pd + 3pd) = ( — 2 + 3n)Gpd- Equation (|T7|) shows 
that Gpd ~ H 2 and Gpd grows with time. If the condition (|22[) is not satisfied, H 
and Gpd diverge in finite time. Furthermore, as pd grows with the scale factor, Gpd 
certainly increases compared to G. For gravitationally bound systems, the GPE 
contribution of the order ~ G(pd + 3pd)R 3 (where R denotes the characteristic 
spatial scale of the bound system) overwhelms the "mass" contribution ~ GM 
{M denotes the mass of the bound system). Gravitationally bound systems fall 
apart in finite time. For the systems bound by electromagnetic or strong forces, 
mere growth of Gpd ensures their unbounding at some finite time before the time at 
which scale factor goes to infinity. Consequently, all bound structures are unbound 
in finite times. The scenario of the "Big Rip" is present in generalized phantom 
energy models as well. 
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Finally, let us make some comments on fundamental aspects of the GPE model. 
As the gravitational constant G(t) is time-dependent, the description of the grav- 
itational sector in the GPE model represents a declination from the Einsteinian 
gravity. One important aspect is whether the scale factor a really describes the 
growth of length scales. One can raise two arguments in favour of the standard 
interpretation of the scale factor a. The first is that no intervention in the geo- 
metrical structure or interpretation of the left-hand side of equation ((7J) has been 
made. The other, more physical one, is that the density of "nonrelativistic" matter 
scales as p m ~ a -3 in our GPE model, equation (j3J) with 7 = 0. Given that no 
interaction (production or annihilation) of the "ordinary" matter component with 
other components is assumed, this fact establishes a as a natural measure of the 
growth of length scales. 

In some theories with the time-dependent effective gravitational constant, such 
as scalar-tensor or nonminimally coupled scalar field theories, one can construct 
many mathematically equivalent theories using conformal transformations. It turns 
out that all these theories are not physically equivalent, i.e. some formulations are 
more physically viable than others (the Einstein frame formulation is more viable 
than the Jordan frame formulation) ^7j- Generally, it might be of interest to 
consider conformally related models of GPE obtained by the transformation of the 
type 9fiv = f(G(t))g^ u , where / is a suitably chosen function. However, the time 
variation of G(t) in our model can be very general and includes possibilities to 
which requirements on the choice of the conformal frame do not necessarily apply. 
Some examples of such a variation are the renormalization group running of G 
[7| |H1 E] or the time variation of G emanating from extra dimensions ^S] . 

In conclusion, in this letter we have considered cosmological models with the 
time-dependent gravitational constant G and dark energy with the supernegative 
equation of state (phantom energy). Phantom energy is generalized in the sense 
that its equation of state does not determine its evolution with the scale factor a, 
i.e. GPE density becomes an independent function of the scale factor. The require- 
ment of general covariance in this setting imposes conditions on the gravitational 
constant G which acquires time dependence. Investigation of future dynamics of 
the generalized phantom energy models with growing generalized phantom energy 
density and the parameter of the equation of state less than —1 exhibits some 
general properties. A large class of models with different evolutions of pd and G, 
but the same equation of state of GPE, have the common law of the evolution of 
the scale factor a in the sufficiently distant future. The time dependence of the 
GPE parameter of the equation of state determines whether the universe evolves 
infinitely in a de Sitter regime or diverges in finite time. One would expect that 
bounds on the variation of G in the past epochs of the evolution of the universe 
would produce the most stringent constraints on the parameters of the GPE model. 
Therefore, it is important to point out that our main results qualitatively do not 
depend on the size of the parameter \n\ or on the intensity of growth of pd (of 
course, within classes of these parameters that satisfy or do not satisfy the con- 
dition (|22|0 . For smaller parameter values and slowlier varying functions pd and 
G, the onset of the general evolution (dependent only on k) will come later. For 
instance, for constant and negative k, but very small \k\, the entire class of GPE 
models leads to the "Big Rip" event, but at very late times. 
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Clearly, the present accelerating phase of the evolution of the universe carries 
the seed of the possibly very dramatic future of our cosmos. Therefore, more 
precise observations of the past variation of and G with time (redshift) will be 
able to unravel the fate of the universe. 
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